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The objective of this effort is to develop an efficient and accurate thermo-fluid 
computational methodology to predict environments for a solid-core, nuclear thermal engine 
thrust chamber. The computational methodology is based on an unstructured-grid, 
pressure-based computational fluid dynamics formulation. A tvro-pronged approach is 
employed in this effort: A detailed thermo-fluid analysis on a multi-channel flow element for 
mid-section corrosion investigation; and a global modeling of the thrust chamber to 
understand the effect of heat transfer on thrust performance. Preliminary results on both 
aspects are presented. 


I. Introduction 

N uclear thermal propulsion (NTP) can carry far larger payloads and reduce travel time of astronauts to Mars than 
is now possible with chemical propulsion. The feasibility of NTP engine systems was established by extensive 
testing in the Rover/NERVA programs and the technical merits of NTP have been identified in numerous studies 
and summarized.’ For the solid-core NTP legacy engines tested during the Rover/NERVA era, heat transfer 
efficiency and degree of hydrogen dissociation are two important performance factors. Yet the impact of hydrogen 
dissociation on engine performance was not quantitatively identified. In addition, although the so-called flow 
element mid-section corrosion^ was thought to be caused by the cracking of coating material and the materials 
development would remove that failure mode, detailed analysis is needed from a fundamental viewpoint to assess 
potential issues and enhancements. The objective of this effort is therefore to develop an efficient and accurate 
multiphysics thermo-fluid computational methodology to predict environments for a solid-core, NTP engine thrust 
chamber, with the emphasis on the effect of hydrogen dissociation on engine performance and investigation of the 
mid-section corrosion mode. 

The current task scope is to develop and implement conjugate heat transfer and porosity modules based on an 
existing Unstructured-grid Navier-Stokes Internal-external computational fluid dynamics Code (UNIC), and to 
perform three-dimensional, multi-physics (real fluid, chemically reacting, turbulent flow, and conjugate heat 
transfer) thrust performance and heat transfer analyses for a solid-core, nuclear thermal engine thrust chamber, to 
provide thermal, fluid, and hydrogen environments. A two-pronged approach is employed in this effort: A detailed 
conjugate heat transfer analysis on a multi-channel flow element, and a global modeling of the thruster chamber 
with a porosity modeling technique. It. is anticipated that the detailed analysis on a single flow element provides 
detailed fluid, thermal, and hydrogen environments for failure investigations, while the global thrust chamber 
analysis helps the understanding of the effect of heat transfer on thrust performance. These modeling activities will 
be validated as much as possible by testing performed by other related tasks. 
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II. Computational Methodology 

A. Computational Fluid Dynamics 

The CFD methodology is based on a multi-dimensional, finite-volume, viscous, chemically reacting, 
unstructured grid, and pressure-based formulation. Time- varying transport equations of continuity, species 
continuity , momentum, total enthalpy, turbulent kinetic energy, and turbulent kinetic energy dissipation were solved 
using a time-marching sub-iteration scheme and are written as: 
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A predictor and corrector solution algorithm was employed to provide coupling of the governing equations. A 
second-order central-difference scheme was employed to discretize the diffusion fluxes and source terms. For the 
convective terms, a second-order upwind total variation diminishing difference scheme was used. To enhance the 
temporal accuracy, a second-order backward difference scheme was employed to discretize the temporal terms. 
Details of the numerical algorithm can be found in Refs 3-7. 

An extended k-e turbulence model® was used to describe the turbulence. A modified wall function approach was 
employed to provide wall boundary layer solutions that are less sensitive to the near-wall grid spacing. 
Consequently, the model has combined the advantages of both the integrated-to-the-wall approach and the 
conventional law-of-the-wall approach by incorporating a complete velocity profile and a universal temperature 
profile®. A 2-species, 3-reaction detailed mechanism’ was used to describe the hydrogen dissociation and 
recombination chemical kinetics. 


B. Computational Conjugate Heat Transfer in Solids 

The solid heat conduction equation is solved with the gas-side heat flux distributions as its boundary conditions. 
The solid heat conduction equation can be written as: 
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where and represent source terms from volumetric and boundary contributions, respectively, ffand C denote 
the thermal conductivity and capacity of the solid material, respectively. The temperature at the fluid-solid interface 
is obtained by enforcing the heat flux conservation condition. 

C. Porosity Modeling 

For the two-temperature porosity model, separate thermal conductivities for the flow and the solid parts are used. 
The heat transfer between the flow and solid is modeled by using the empirical correlation of heat transfer 
coefficient for circular pipes as a function of flow Reynolds number. Ergun’s equation is used as a point of departure 
for the drag loss. An empirical multiplier for both the heat transfer and drag loss will be determined by comparing 
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solutions of flow passing through a porous flow element with those of a 19-chaimel flow element using detailed 
conjugate heat transfer modeling. The Continuity, Navier-Stokes and Energy (Total Enthalpy) Equations, can be 
written in a Cartesian tensor form with a porosity factor, a, included: 
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For the solid heat conduction, 
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where and Qa represent the drag loss and heat transfer source due to the porous material. Using the Ergun’s 
equation for the drag loss. 
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For the heat exchange source term, 

Q^=fA^H[T,-T^) (13) 


where H and A„ represent the heat transfer coefficient and wetted surface area between the gas and solid, 
respectively. Subscripts s and g denote the solid and gas domain, respectively. The remaining parameter on the 
right-hand-side of Eq. (13),/, is an empirical parameter to be determined by comparing solution of flow passing 
through a 19-channel fuel element using detailed conjugate heat transfer modeling and that of flow passing through 
a homogeneous fuel element with the effect of the 19-channel described by the porosity model. 


III. Computational Grid Generation 

Hybrid computational grids were generated using a software package GREDGEN. A series of grid verification 
studies were performed to determine the current grid size. It is also found during the grid study that the computed 
gas-solid interface temperature is most stable when structured-grid layers are present on both sides of the interface. 


IV. Boundary and Inlet Conditions 

No-slip condition was applied to the solid walls. Fixed mass flow rate boundary condition was used at the inlet, 
and mass conservation boundary condition was used at the exit. The inlet conditions were obtained from system 
simulation for the Small Engine. 


V. Results and Discussion 
A. Preliminary power distribution for a single flow element 
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It is known that the power distribution inside a reactor normally resembles a Cosine curve with its peak at the 
center and reduces to a small value at the boundaries due to escaping neutrons. To model this effect numerically for 
a single, chamber-sized flow element, a section is implemented in the porous volume identification subroutine to 
define the boundary distances from the axis of peak power generation. This set of distances is then used to construct 
the Cosine distribution function. This method is good for fuel assembly regions with constant cross-sectional shape. 
An off-centered peak power input can be assigned. 

Figures 1 and 2 show the predicted chamber temperature contours for a centered and an off-centered peak power 
distribution in a single porous flow element. The peak chamber temperatures for the two cases are about the same, 
but the location of the peak chamber temperature of the off-centered power distribution solution is shifted towards y 
= -2 cm location. Due to this shift in temperature profile, the density and velocity distributions also become biased. 
This in turn causes biased drag loss and pressure distributions in the porous region. As a result, side force is created. 
It is postulated that biased power distributions can end up with substantial side forces, which could be a contributing 
factor in flow element failures. 
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Fig. 1 Predicted chamber temperature contours with 
centered peak power in the porous flow element 
region from 3 to 10 cm. 
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Fig. 2 Predicted chamber temperature contours 
with peak power off-centered (along y = -2 cm 
and z = 0 cm axis). 


B. Initial 3-D grids for a detailed 19-channel flow element 

Figure 3 shows the surface grid representation for the web of a 19-channel flow element. Figure 4 shows a 
detailed grid representation of the web and hot flow channels. Flow solutions for a series of fixed outer wall 
temperatures will be performed in order to validate the heat transfer equation for the porosity model. Flow solutions 
for different channel diameters at 2.03, 2.29 and 2.54 mm will also be performed to anchor the drag law for the 
porosity model. 
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Fig. 3 Surface grid representation for the web 
of a 19-channel flow element. 



Fig. 4 Detailed grid representation of a quarter 
of a fuel element. 


VI. Conclusions 

An efficient and accurate multiphysics thermo-fluid computational methodology is being developed to predict 
environments for a solid-core, NTP engine thrust chamber, to study the effect of hydrogen dissociation on engine 
performance and investigation of the mid-section corrosion mode. A two-pronged approach is employed in this 
effort: A detailed conjugate heat transfer analysis on a multi-channel flow element, and a global modeling of the 
thruster chamber with a porosity modeling technique. It is anticipated that the detailed analysis on a single flow 
element provides detailed fluid, thermal, and hydrogen environments for better understanding of the mid-section 
corrosion phenomenon, while the global thrust chamber analysis helps the understanding of the effects of heat 
transfer and hydrogen dissociation on thrust performance. 
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